RELAXATION TIME OF ANISOTROPIC SIMPLE EXCLUSION 
PROCESSES AND QUANTUM HEISENBERG MODELS 
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Abstract. Motivated by an exact mapping between anisotropic half integer spin 
quantum Heisenberg models and asymmetric diffusions on the lattice, we consider an 
anisotropic simple exclusion process with N particles in a rectangle of Z 2 . Every particle 
at row h tries to jump to an arbitrary empty site at row h ± 1 with rate q , where 
q 6 (0, 1) is a measure of the drift driving the particles towards the bottom of the rectan- 
gle. We prove that the spectral gap of the generator is uniformly positive in N and in the 
size of the rectangle. The proof is inspired by a recent interesting technique envisaged 
by E. Carlen, M.C. Carvalho and M. Loss to analyze the Kac model for the non linear 
Boltzmann equation. We then apply the result to prove precise upper and lower bounds 
on the energy gap for the spin-S, S £ |N, XXZ chain and for the 111 interface of the 
spin-S XXZ Heisenberg model, thus generalizing previous results valid only for spin 1 
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1. Introduction 

Some years ago it was discovered by Alcaraz [2] that a class of asymmetric reversible 
simple exclusion processes on Z rf related to models of diffusion limited chemical reactions, 
are unitarily equivalent to certain anisotropic quantum Heisenberg Hamiltonian, known 
as XXZ models, that have received in recent years increasing attention in connection with 
the analysis of quantum domain walls (see [1], [11], [14] and references therein). Such 
an equivalence implies that the spectrum of (minus) the Markov generator of the process 
coincides with the spectrum of the quantum Hamiltonian. In particular the energy gap 
above the quantum ground state, a key quantity in the theory of quantum spin systems, 
becomes identical to the spectral gap of the process and a variety of probabilistic techniques 
come into play in order to obtain meaningful estimates. Such an observation was exploited 
recently in [7] to prove sharp bounds on the energy of low lying excitations above a ground 
state describing a 111 interface for a spin \ XXZ model. The extension of such results to 
higher half integer spin S 6 requires additional finer analysis and led us to consider 
the following model. 

1.1. Setup. Given two natural numbers L,H we consider the rectangle 

A = {(i,/i)GZ 2 : i = l,...,Landh = l,...,H} (1.1) 

For each i, Aj stands for the stick at i given by A, = {(i, h) : h = 1, ... ,H}. At each 
x € A we have a variable ot x £ {0, 1}: we say that site x is occupied (by a particle) if 
ol x = 1 and empty otherwise. The set of configurations {0, 1} A is denoted by f2 and it is 
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naturally decomposed in single stick configurations: a G fi will be written often in the 
form a = (771, ... , tjl) with rji £ {0, 1}^ denoting the restriction of a to the stick Aj. 
Given a parameter q £ (0, 1) we define the product probability measure fx on fi. 

m(/) = e , M«)=nnfri ( L2 ) 

a€S] i=l h=l ^ 

where / is a generic function / : $7 — ► R. According to /U particles prefer to live on the 
region of small h, i.e. the bottom of the box A if we interpret h as a vertical coordinate. 
We define nj as the number of particles in the stick Aj: rii(a) = rii(r]i) = J2h=i a (i,h) an d 
consider the conditional probability measure 

L 

v = y N = fx( . \y^m = n) (1.3) 

i=l 

The variance of a function / w.r.t. u will be written as usual in one of the following ways 

Var(/) = K/,/) = K(/"K/)) 2 )- 

1.2. The process and main result. The asymmetric diffusion that will be analyzed 
in the sequel can be described as follows. Every particle at row h tries to jump to an 
arbitrary empty site at row h + 1 with rate q and to an empty site at row h — 1 with rate 
1/q. The Markov generator is defined by the operator 

j L L H-l 

= 1 E E E c (i,h);(j.h+l)( a ) V (<.'»);(j.'H-l)/( Q! ) ( L4 ) 
i=l j=i /i=i 

where we use the notation 

V(^ );0>+1) /(a) = /( a (^);0>+D) _ /(a) , (1.5) 

q,(»,/i);0',/i+i) denoting the configuration in which the values of a at (i, h) and (j, h+ 1) have 
been interchanged while the rest is kept unchanged. The rates t h);(j,h+i) are given by 

c ( i 1 h);( J - 1 h+i)(«)=9 a(il ' ,) " a(j ' I ' ,+1) - (1-6) 

Simple computations show that C is self adjoint in L 2 (u), the associated Dirichlet-form 
being 

L L 

1=1 i=l 

1 2 

h=l 

Our main result then says that decay to equilibrium for the dynamics defined by (1.4) 
occurs exponentially fast in the L 2 (u) norm, uniformly in L,H and N. As a corollary 
we shall obtain an interesting estimate on the energy gap for a class of quantum XXZ 
Hamiltonian, see Theorems 4.1 and 5.1 below, that extends in particular previous results 
in [7] and [10]. 
More precisely let 

7 (L,#) = sup sup (1.8) 
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where the number of particles N in sup^y, using the particle-hole symmetry, is assumed 
to range from 1 to 

Theorem 1.1. For every q G (0, 1) there exists C < oo such that 

sup7(L,iT) C . 

L,H 



Remark 1.2. Recently we learned [3] that for L = 1 and N = If the mixing time (namely 
the smallest time such that sup^ v (\e tc (a, •) — e tc (a, -)|) < ^) grows like H. Remarkably 
in the same setting the logarithmic Sobolev constant grows at least like H 2 . 

1.3. Applications. Some of the applications of Theorem 1.1 we have in mind, particularly 
those to quantum Heisenberg models, are linked to the analysis of the restriction of the 
above defined process to the horizontal sums of the basic variables given by 

L 

Uh = ^2 ®(i,h) , h = l,...,H 
i=i 

We will show below that the evolution of the new variables {o->h}h=i ^ s still Markovian and 
that it can be interpreted as describing the fluctuations of a non-negative profile subject 
to a fixed area constraint. 

Let Vl denote the set of permutations of {1, . . . , L}. Given it G Vl we write a n,h for the 
configuration 

p ; _K,fc') h '* h , i = i,..., L (1.9) 

The subspace S of horizontally symmetric functions is defined by 

S = {f G L 2 (u) : f{a) = f(a*> h ) , Vtt G Vl, = 1, . . . , H} . (1.10) 

and it clearly consists of functions which only depend on the horizontal sums. 
Given / G S we write f(i~o) = f(a). In this way we identify S with the space L 2 (r2,z>), 
(l = {0, 1, . . . , L + 1} H and z> the marginal of v on horizontal sums u = {u^}. The 
probability i>(w) of a single u> G & compatible with the global constraint Ylh^h = N is 
easily computed to be 

ojGfi: ^— 1 

Observe that S is an invariant subspace for the generator C, i.e. CS C S. In fact a simple 
computation shows that for every / G <S, Cf(a) can be written as 

JZ-l 

= i E { ™+,*M [/V'") - /H] + t/V ,fc ) - /H] } (i-i2) 
h=i 

w+,h ■= q' 1 (L - u h )u h+1 , w-, h :=q(L- u h+1 )u h , 

ujh' h! 7^ h, h + 1 

w^,'' 1 := < a;^ ± 1 h! = h 

w W Tl h' = h+l 
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This defines a Markov generator C which is symmetric in L 2 (Q,u). The corresponding 
process (the restriction to {uJh} of the original anisotropic exclusion dynamics) can be 
interpreted as describing fluctuations of a non-negative profile u> := {^h\h=i subject to a 
fixed area constraint (^^oj/, = constant). In view of the anisotropy the profile is rather 
sharply localized: letting p = N/L we see that ~ L for heights h below p and ~ 
above p with high probability. By Theorem 1.1 relaxation to equilibrium in L 2 (&,£>) is 
exponentially fast uniformly in p. 

In the case L = 2 the previous analysis admits another interesting interpretation as a 
model for diffusion limited chemical reactions, see [2, 5, 6] and references therein. Namely 
describe the state u>h = 2 as the presence at h of a particle of type A, ujh = as a particle 
of type B and u>h = 1 as the absence of particles (inert). If tia, ub denote the size of 
the two populations we see that the difference tia — n# is conserved and this system can 
be studied as a model for asymmetric diffusion with creation and annihilation of the two 
species. Particles of type A have a constant drift towards the bottom ( "small h" region) 
while particles of type B have the same drift towards the top ("large h" region). They 
perform asymmetric simple exclusion with respect to the inert sites but when they meet 
(i.e. when they become nearest neighbors) they can produce the annihilation reaction 
A + B — > inert. The reverse reaction inert — ► A + B restores steady state fluctuations 
given by the canonical measure. 

While Theorem 1.1 implies immediately L 2 -exponential ergodicity for the above process, 
a direct proof of the result for the two-particle model seemed difficult to us. 

1.4. Some ideas for the proof of Theorem 1.1. We conclude this introductory section 
with some comments on the main ideas behind the proof of Theorem 1.1. Our main source 
of inspiration has been a recent work by E. Carlen, M.C. Carvalho and M. Loss on the rate 
of approach to equilibrium for the Kac model of the non linear Boltzmann equation [8]. 
Like in other approaches to bound the spectral gap for large reversible Markov chains, the 
first idea is to recursively bound 7(L, H) in terms 7(1, H), the latter being finite uniformly 
in H because of Theorem 4.3 in [7]. The starting point, as e.g. in the martingale approach 
of H.-T. Yau [15], is a decomposition of the variance of an arbitrary function / as 

L L 

Var(/) = -Y, f(Var(/ I ?k)) + j E Var (K7 I ^)) • 
fc=i fc=i 

where J~k denote the a— algebra generated by the stick-variables k = 1, . . . , ,L. It 
is easy to check (see section 3 below) that the first term can be bounded in terms of 
j(L — 1, H) x V(f, f). The main new idea comes in the analysis of the second term and 
consists in introducing the stochastic symmetric operator 

1 L 

fc=i 

and observing that for any mean zero function / the following identity holds true: 

1 L 

jY i -Vax{i,tf\r k ))=v{fPf). 

fc=l 

Thus 

1 L 

Var (/) - 1 E Var ("(/ I ^)) = Var(/) - u(fPf) = v(f(l - P)f) . 
k=i 
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so that one is left with the problem of establishing an estimate from below on the spectral 
gap of P which is sharp enough to allow a successful iteration in L for j(L,H). The 
key point now is that, because of the particular form of P and of the symmetry of the 
measure the estimate of the spectral gap of P boils down to the estimate from below 
of the spectral gap of a particular one dimensional random walk that can be described 
as follows. Let n± be the minimum and maximum number of particles allowed in a 
single stick, say the first one. Then the state space for the random walk is the interval 
[n_, n_ + 1, . . . , n+] and the transition kernel q(n — > m) is given by v(n\ = m \ n<i = n). 
It is easy to check that such a process is ergodic iff L > 3. The study of its relaxation 
time represents in some sense the technical core of the paper and it requires a rather fine 
analysis based on equivalence of ensembles type of results. 
The rest of the paper is organized as follows. 

(i) In the next section we define and analyze the one dimensional random walk mentioned 
above. 

(ii) In section 3 we prove Theorem 1.1. 

(iii) In sections 4,5 we discuss the main applications of Theorem 1.1 to quantum Heisen- 
berg XXZ models. 

Acknowledgments. Part of this work was done at the Institute H. Poincare during 
the special semester on "Hydrodynamic limits". We would like to thank the organizers 
F. Golse and S. Olla for their kind invitation and the stimulating scientific atmosphere 
there. We are also grateful to T. Koma, B. Nachtergaele and S. Starr for informing us 
about their results prior to publication and to I. Benjamini, N. Berger, C. Hoffman and 
E. Mossel for an interesting discussion concerning the mixing time of asymmetric simple 
exclusion. Finally we would like to thank F. Cesi for a very interesting and enlightening 
discussion which helped us to clarify a tricky point in our argument. 
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2. Spectral gap of the long-jump random walk in a single stick 

In this section we are going to study a one-dimensional process which plays a key role 
in the recursive proof of Theorem 1.1. 

Let uq be the marginal of v, the canonical measure defined in (1.3), on a single stick 
configuration. If 7Tj : {0, 1} A — ► {0, l} Ai denotes the canonical projection onto single 
stick configurations {jtiOi = rji) we may write uq = v o ir^ 1 . Let EI denote the space 
L 2 ({0, l} Al , vq) and denote by (■, •) the corresponding scalar product: 

{(p,tp) = v((ip o7ri)(^o7ri)) , p,ipeW. 

The ex-algebra generated by the single stick variable rji is denoted by Ti. The operator 
K : H -> H is defined by 

Klf = v(p O 7T 2 | J^l) , 

or equivalently by the bilinear form 

{if, Kip) = v(((f O 7Tl)(V> ° 7T 2 )) . 

JT is a stochastic, symmetric linear operator on H. The number of particles in a stick 
is denoted by n and we call p = N/L its average value. The centered variable n — p is 
denoted by n. Observe that n is an eigenfunction of K with eigenvalue — 1/(L — 1): 

Kh = -- -n. (2.1) 

lj — l 

Our main result in this section states that apart from the values 1 and —\/{L — 1) the 
spectrum of K is concentrated around within an interval of radius 0(L~ 1-5 ) for some 
5 > 0, uniformly in TV, H. 

Theorem 2.1. There exist constants 5 > 0, k < oo and Lq < oo independent of N,H 
such that if L ^ Lq 

\( i p,Kp)\^kL- 1 - s ( i p,p) (2.2) 
/or aZZ <£> G H with vo(<p) = and (99, n) = 0. 

Proof. Let Ho denote the subspace of H of functions which only depend on the number 
of particles. Let also E denote the orthogonal projection onto Ho, i.e. if J-q denotes the 
o"-algebra generated by the random variable n we have 

Ep = z^o I -^b) • 

A simple computation now shows that K commutes with E and 

Kp = KEp , p e H . 

These observations prove that (tp, Kp) = {Ep, KEp) and in order to prove the claim (2.2) 
we can restrict to p £ Mq. 

Let us introduce some handy notations. We write v(n) for the probability that in one 
given stick there are n particles, v(nC\m) for the probability that in two given sticks there 
are n and m particles respectively, and v{n \ m) for the probability that in one given stick 
there are n particles conditioned to the event that in a different given stick there are m 
particles. With these symbols we have 

Kp(n) = u ( m I n)p{m) p G Ho . 

m 
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The operator K then describes a random walk on the integers with transition probabilities 
Prob {n — > m) = v(m\n). The desired estimate (2.2) can be written as 

v(m)v(n)ip(m)Q(m, re)y(n) ^ k L^ 1 ^ 6 {f,<p} , (2-3) 

«(»,»):=^-l (2.4) 

for all 99 G Ho such that fo(^) = and (<£>, n) = 0. The rest of this section is concerned 
with the proof of (2.3). 

The idea is to split the sum in (2.3) in a region of typical values of m, n where things 
are controlled by a careful expansion and a region of atypical values whose contribution 
is shown to be negligible by tail estimates. Unfortunately the definition of typical and 
atypical values of n, m strongly depends on the value of the particle density p and we will 
be forced to distinguish between two cases, conventionally denoted large density and 
small density case, depending whether p > L~4 or p < L~*. Before entering into the 
details of the proof we first need to establish some preliminary useful bounds. 

2.1. Technical bounds. The grand canonical measure p = p x ^ with density p = N/L 
is the product measure on A (see (1.2)) with 

^)=nn + ^» (2 - 5) 

i=lh=l H 

Here the parameter X(p) £ R, often called the chemical potential, is such that the average 
number of particles in any given stick is equal to p. The variance p(n 2 ) of the number of 
particles in a stick is denoted by a 2 . Simple computations - as in [7], Lemma 3.2 - show 
that for every q G (0, 1), there exists k < 00 such that 

k~ x {p A 1) ^a 2 ^k(pM). 

Similarly p(n 4 ) ^ k(p A 1). As a rule, here and throughout the rest of this section the 
letter k will be used to denote a finite constant whose value may change from line to line. 
What is essential is that it only depends on q and is uniform in all other parameters: L,H, 
and N. 

We introduce the characteristic function 

F(t) = p(e l ^L n ) . 
We shall rely on the following simple estimate. 
Lemma 2.2. For all t G [— ira\^L, ■KO\fT\ 

\Fit)\ ^ e~ ki r (2.6) 

Proof. Writing 

n = ah ' 

h 

where we drop the horizontal stick label, we have 



\F(t)\=l[g h (t), g h (t):=\p(e l ^T ah ) 
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Using the inequality x ^ e^ 2 l ^ 2 , x G (0, 1), we have 
g h (t) ^ exp {\(g h {tf - 1)) 

= exp{ - ^Var M [cos (ta h /aVI)] - ^Var M [ sin (ta h /aVI)} } 
Define a 2 = Var /J (a/ l ) and compute 

Var^[cos (ta.h/a\fL)\ = ct^cos (t/aVZ) — l) 2 
Var M [ sin (ta h /aVL)] = o\ ( sin (t/aVT,)) 2 . 

It follows 

g h (t) ^ exp { - o\{\ - cos (t/a^L))} ^ exp ( - kt 2 a 2 h /a 2 L) , 

where we use the inequality 1 — coss ^ ks 2 , for some k > and all s G [— 7r, it]. The 
lemma now follows from a 2 = Ylh a \- '-' 

It is useful to compare the canonical measure v, given again by (1.3), with the grand canon- 
ical probability p. In particular, it follows from Lemma 2.2 - as shown in [7], Proposition 
3.8 - that there exists k < oo independent of the density p such that 

u(f)^k^(f) (2.7) 

for any function / > depending at most on ^ variables {r/j}. We shall need the following 
estimates on the tails of our distributions. 

Lemma 2.3. There exist constants a > and k < oo depending only on q such that 

v(n | m)v(m \n)<ke' a l(n- P ) 2 +(™-p) 2 ] ( 2 .g) 
v{n\m)v{m\n) < k p n+m (2.9) 

Proof. We start with the proof of (2.8). Given n G N we denote by X(n) the chemical 
potential such p x ( n \ni) = n. We write simply A for the chemical potential \(p). We have 

v{n) < kp x (n) = ke^ (n) p x(n \n) < ke~ v ^ (2.10) 

where 

/ £A(n)\ 

V\(n) = c (A(n) — A)n — log I — tft - I . c:=-logg, 
and Z x denotes the partition function 

^ (l + q 2{h ~ X) ) ■ 
h=l 

Now V\(p) = and we write 



A standard computation gives that j^V\(s) = c(A(s) — A) and 

^ Vx{s)=c i x{s) = ^m^) >a 
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for some a > independent of s. Therefore 

V x (n) >a(n-pf. 
By (2.10) we have the tail estimate 

v(n) < ke- a ^" p ^ . (2.11) 

We now apply the above bound to each term v(n \ m). In this case the starting chemical 
potential A corresponds to a density 

p' = tEt so that (2- 11 ) y ields 

„ (_ N-m \ 2 

u(n\m) < ke~ a y n ~-^> 
Our original claim (2.8) now follows from the uniform (in N, L) convexity of the function 

N-y\ 2 ( N-x^ 2 



G(x,y) := [x - — ) + \y ^ — 1 
around the unique minimum x = p, y = p. 
To prove (2.9) we use again the bound (2.7) and estimate 

u(n\m) < kp [ -^>(n) < k ( — — -J < hp n . 



□ 



We are in a position to complete the proof of Theorem 2.1. As anticipated we will need 
to distinguish between "large" and "small" values of the density p. 

3 

2.2. The large density case. Here we assume p ^ L~i. Define the set 

B = {(n,m) : |n| + \m\ ^ BlogL}, (2.12) 
with B < oo a constant to be fixed later. We rewrite the sum in (2.3) as 

v{n)v{m)(p{n)Q{m,n)Lp{m) + v(n)i'(m)(p(n)Q(m,n)(p(m) . (2.13) 

(n,m)eB (n,m)£B 

The second term in (2.13) is estimated with the help of Lemma 2.3. To see this use 
Schwarz' inequality to write 

v(m)v(n)(p(m)Q(m,n)(p(n) 

(n,m)£B 

\2 



< (ip,(p) ^ v(m)v(n)Q(m,n) 

(n,m)£B 



<(<p,<p)[ E " (n)i/(m)2 [1 + ^ 



v{n | m) 



2 



(n,m)<£B V ; 



v(n)v{m) — ^ — = v{n \ ra)v(m \ n) 



Since 

v(n | m) 2 

an application of (2.8) and (2.11) gives 

v(n)v(m)<fi{n)Q(m,n)(p(m) ^ k L~ 2 (<f,(f) , (2-14) 

n,m^B 
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provided B is sufficiently large. 

The key step in the proof of Theorem 2.1 in the case of large density will be the following 
expansion. 

Lemma 2.4. For all q G (0, 1), B < oo, there exist constants k < oo, ( > such that if 

_ 3 

p ^ L * then 



nm 

Q(m,n) = — + K(n,m) (2.15) 



with the remainder 1Z satisfying 



1 1/ 2 

v(n)v(m)\K{n,m)\ 2 < kL' 1 ^ . (2.16) 



n,m 

\n\ + \m\ ^ BlogL 



Proof. We write 



u{n) = ^ f dtF{t) L ~ l e^ n 

v(nf\m) = ■= / dtF[t) e 

2tto-\/ Lfi(N) J 

where all the integrals are over the interval [— iras/Z, jrcry/L]. Since 

2 7 rav / I/i(iV) = J dtF(t) L 



we can write 



with 



and 



v(n PI m) — v(ri)v(m) NUM 
Q(m,n) = u(n)u(m) = DEN 



NUM := JdtF(t) L - 2 e l ^ [n+fh] J dt! F(t') L 
-J dtF{t) L ~ l e l ^L n j dt! Fitlf^e 1 ^ 

DEN := J dtF(t) L - 1 e i ^Z ii J dt! Fitlf^e 1 ^ 



Notice that, because of the Gaussian upper bound of Lemma 2.2 we have |-F(t)| L ^ e" at 
and only the region \t\ < klogL (for some large but fixed k) will have to be taken care 
of. In order to be precise about the nature of the various error terms, in what follows we 
will denote by e(L) a generic term which, upon multiplication by L, still goes to zero (as 
L — > oo) as an inverse power of L uniformly in the range of \t\ ^ ATogL. We first observe 
that F(t) = 1 — + e(L). Indeed by expanding F around t = the third order error 
term is bounded from above by 

k jMw~ 2 ^ ) - k ^ - £{L) 

where we use the bound 

/J.(\n\ 3 ) < fccr/i(n 4 )5 ^ ko 2 
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together with a 2 ^ k(p A 1) and p ^ L 4 . This implies F(t) 1 = 1 + + e(L) and 
F(t) -2 = 1 + £ + e(L). Then if we write e i ^T fi = 1 + <J n (t) we have 

NUM = y dtF(t) L 5 n (t)5 m (t) J dt'F{t') L - j dtF(t) L 5 n (t) J dt'F(t') L 5 m (t') + e(L) . 

If we modify further and define 



S n = 5 n -i 



n . 



I(L) = J dtF{t) L t 



we then have, 



NUM 



nm 



a 2 L 



+ i 



dtt 2 F{t) L j dt' F{t') L - /(L) 2 ) 
+ J dtF(t) L S n (t)S m (t) I dt'F{t') L - [ dtF{t) L 8 n {t) [ dt' F{t') L 5 m {t') 

— I ill wi i\ L 1 A ( 1 \ I ,11' rn'\ L 4 ; " / ,// n ; J- 1 :< 1 r\ l,n'^(i') L 
u y/~L 



dtF{t) L t5 n (t) J dt' F(t') L + i-^= J dtF{t) L t5 m {t) J dt' Fit') 1 



m 



J dtF(t) L 5 n (t) J dt'F(t') L t' 



n 



Now observe that \F\V) = e(L) and 



dt 



F{t) L 5 m {t) J dt'F(t') L t' + e{L). 



i j dt t 2 F(t) L J dt' F(t') L = ^ + e{L) . 



Using also \6 n (t)\ < k rj -£- it follows that 



<T 2 L 



nm 



NUM = -2vr — - + R(n, m) + e(L) 



with 

\R(n,m)\ < k 
What is crucial for us is that 



\n\ \m\ 



<7" 



e(L) + 



n 2 m 2 \n\ m 2 + \m\ n 2 



+ 



(a 2 L) 2 {a 2 Lf/ 2 



1/2 



u(n)u(m)\R(n,m)\ 2 < e(L) 



(2.17) 



\n\+m\ ^ B log L 



provided that p > L~3. The above estimate actually holds without the restriction |n| + 
\m\ ^ BlogL as it is easily seen using z^(n 4 ) ^ fc^(n 4 ) ^ ho 2 and the bound (pAl) ^ ho 2 . 
On the other hand if we repeat the reasoning for the denominator DEN we obtain 

DEN = 27r + R(n,m) 

with a remainder R(n, m) satisfying 

sup \R(n,m)\ ^ kL~^ 



n,m: 
\n\+\fh\ < BlogL 



(2.18) 



for some ( > whenever p ^ L 4. In conclusion Q(m,n) has been written as in (2.15) 
and (2.17)-(2.18) imply (2.16). □ 
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We are now able to finish the proof of Theorem 2.1 in the case of large densities. Observe 
that 

u{n)u{m)(p{n)(p(m)nm={if,n) 2 — v(n)v(m)<p{n)ip{m)nrh . 

(n,m)eB (n,m)£B 

By assumption (ip, n) = 0. Then 

v(n)v(m)Lp{n)ip{m)nfh 

(n,m)£B 

^ |^ v(n)v(m)n 2 fh 2 2 {<p,<p) ^ kL~ 2 {<p,<p) 

where the last estimate can be easily obtained from (2.11). By Lemma 2.4 we then have 
u{n)u{m)if(n)Q{n, m)ip(m) | 

(n,m)GB 

sC (<p,cp){kL- 2 + [ ^ ^{n)u{m)\TZ{n,m)\ 2 ^ /2 ^ ^ k L -1 ^ (p )¥ >) 

3 

Together with (2.14) these bounds imply Theorem 2.1 when p ^ L~~z. 

3 

2.3. The case of small density. When p < L~i the strategy for the proof of (2.3) 
has to be slightly modified since various previous technical estimates are no longer valid. 
We can however take advantage of the thinner tails of the distribution of the number of 
particles. In this respect we observe that 

i 

v(n)i<(m)Q(n,m) 2 2 < kL' 1 ^ (2.19) 

n,m 

n+m ^ 2 , nm^l 

with C = 1/8. Indeed, thanks to (2.9) 

i 

^2 v(n)v(m)Q{n,m) 2 2 < k pi k L _1_c . 

n,m 
n+m>3 



On the other hand we can examine explicitly the case n = 2, m = 0. We have 
Q(2 Q)- 1 1 = 2} n {m > 1» 

Using the bounds u(0) ^ 1 — kp, v(n) ^ kp n and v({n = 2} n {m > 1}) ^ kp 3 one has 

K2)K0)Q(2,0) 2 < k P 4 - 

We have proved (2.19). 

At this point we may proceed as in (2.13) with the choice 

B := {n, m : n < 1, m < 1} . 

The second term in (2.13) is controlled by (2.19). The first term is given by 

y^ ip(n)ip(m) [ v(n n m) — v(n)v(m) ] . 

n<l m<l 
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We are going to study the eigenvalues of the symmetric 2x2 matrix 

At(n, m) := [ v(n n m) — v(n)v(m) ] , n, m G {0, 1} . 

For n = 0, 1 . . . N, let M n be the number of sticks with exactly n particles. The identities 
L = En M n and iV = E„ imply 

M = L - N + J^(n - l)M n , Mi = iV - ^ nM n . 

n>2 n>2 

Denoting by rii the number of particles in the stick Aj, we have the estimates 
Var^(Mo) < War„ (^(m - l)I {n . >2} ^ < kLp 2 

Var„(Mi) < War„ ^>;I {n . ^ < kLp 2 . 

We compute 

1 L 

i/({m = 0} n {n 2 = o}) = — - "({ni = 0} n {m = o}) 



i=2 



)) 



z,(l {ni=0} (M - 1)) = — _ i/(M (M - 1 

Var(Mo) + — - - ^(M ) 2 - — - - u{M ) . 



L(L - 1) u/ L(L- 1) v VJ L(L — 1 

It follows that At (0,0) = i/({ni = 0} n {n 2 = 0}) - v(l)u(0) can be written 



Similarly 



M(0, 1) = + Cov,(M , MO . 



Introducing the matrix 

the above computations show that 



-1 1 
1 -H ' 



At = y^-A + At 
i/ — i 

with a symmetric matrix At such that |At(n, m)\ ^ 0(p 2 /L), n,m G {0,1}. Then we 
have 

| ^ ^(nMm)At(n,m)| ^^[^0)-^l)] 2 + ^[^(l) 2 + ^(0) 2 ] (2.20) 

n < 1 
m ^ 1 
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Observe that by the orthogonality (if, n) = and Schwarz inequality 
v{n)v{m)(p(n)(p(m)nfh 



n < l 

m < 1 



i 

^ (<p,(p) | v(n)v(m)n 2 fh 2 2 ^ k p2 {(p,tp) 



n+m ^ 2 



On the other hand 

!/(n)i/(m)p(n)^(ffl)iim = p 2 [^(0) - </?(l)] 2 + 0{ P z )[ip{lf + ^(O) 2 ] . 

n ^ 1 
m ^ 1 

Collecting these estimates we arrive at 

[p(0) - ^(l)] 2 < fcp-s (99, tp) + kp[ V (l) 2 + 99(0) 2 ] 
Using also the bound 

Mlf + ^O) 2 ] Op" 1 <¥>,¥>> 
and going back to (2.20) we finally obtain 



| ^ <p(n)<p(m)A1(n,m) < k^j- (<p,(p) ^ kL 1 ^ (<p,<p) 

n ^ 1 
m ^ 1 

This ends the proof of Theorem 2.1 
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3. Proof of Theorem 1.1 

Recall the definition of the stochastic operator K introduced at the beginning of the 
previous section. We set 

w{L,H) =sup[gap(!-i<Q]~ 1 (3.1) 



N 



where 



Define also 



8ap(I-J0= mf i^i^M. (3.2) 

^GH: {if, if) 



7(L,#) = sup sup (3.3) 

N feL*(y) T>(f,f) 



where we introduce the modified Dirichlet form 

i j^i 

Clearly V(f,f) ^ V(f,f) for all / and j(L, H) < j(L,H). Note also that V(f,f) is 
still ergodic: T>(f, f) = implies / = const, for all L ^ 3. More precisely, using e.g. the 
method of Lemma 2.6 in [7] one can easily prove that there exists k = k{q) < oo such that 
for every L ^ 3 one has V(f, /) ^ kV(f, f) for all / € L 2 (v). The key step in the proof of 
Theorem 1.1 is represented by the following proposition whose proof is postponed to the 
end of the section. 

Proposition 3.1. There exists Lq < oo such that for any L ^ Lq and any H ^ 1 

y(L,H)^[lVw(L,H)]%L-l,H) (3.5) 



To complete the proof of Theorem 1.1 we need an estimate on w(L, H). In view of Theorem 
2.1 we know that for every (p G H, such that vo(<p) = we have 

{<p, (1 - K)<p) >{l-k L- 1 - 6 ) (<p, if) (3.6) 

for L ^ Lq with uniform constants 5 > and Lq, k < oo. It is then immediate to deduce 
w(L, H) ^ 1 + kL- 1 ' 5 and therefore 

oo 

l\w(L,H)^C, (3.7) 

L=L 

with some uniform constant C < oo. Then Theorem 1.1 follows from Proposition 3.1 and 
the bound 

sup7(Lo, H) < oo , Lo ^ 3 . 

H 

The latter is easily deduced e.g. from results in [7]. 
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3.1. Proof of Proposition 3.1. Let T k denote the a— algebra generated by the stick- 
variables 7]k, k = 1, . . . , , L. For any / one has the decomposition 

1 L 

Var (/) = 1 E {"( Var (/ I ^)) + Var (K/ I } • (3.8) 



k=i 

We first establish the estimate 

L 



j J2 ^(Vai(/ | F k )) < 7 (£ - 1, #)£(/, /) • (3.9) 

k=l 

By definition of 7 for every /c one has 

^(Var(/ I F k )) ^(i-l^^EE • 
Summing over and using 

E E E °M = ( L - 2 ) E E ^(/) = L ( L - 2 w» /) 

we obtain (3.9). 

We turn to an estimate on the second term in (3.8). Consider the non-negative stochastic 
operator P : L 2 (v) — > L 2 (u) defined by 

1 L 

p / = iE^/i^)- ( 3 - 10 ) 

k=l 

We may assume without loss that u(f) = 0. Then we have the identity 

1 L 

-Y,Vax{vtf\r k ))=v{fPf). (3.11) 

k=l 

From (3.8), (3.9) and (3.11) we obtain 

i/(/(I - P)f) < ^| 7(£ " 1, #)£(/, /) • (3.12) 

In order to estimate from below the left hand side in (3.12) we are going to prove the 
bound 

gap(I-P) ^ ^| [lAgap(I-K)] . (3.13) 

Here gap(l — P) stands for the smallest nonzero eigenvalue of 1 — P. Note that (3.13) 
and (3.12) immediately imply the proposition. 

Take / G L 2 (y) such that v{f) = and 

Pf = Xf, A>0. 

Then / is of the form 

/fa) = (3-14) 



with ^ei and we have the identity 

= L 



^^2Mne) = j^2^2^(fk(m)\ve) ■ (3.15) 

1 k 
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Define now the function $j £l: 

$ / = Yl w ■ 

i 

Taking conditional expectation with respect to Tj in (3.15) a simple computation yields 

\K$ f + A(l - K)<pj = ^-^-K 2 <5> f + \k{2 - K)Q f + y(I - Kfifj . 
Summing over j and factorizing we obtain 

$/ = (3.16) 

The above identity says that if (K + x~t)^/ / then ^ := i s i n the spectrum of K. 
In this case then 

l-\ = ^l{l- li ) > ^-lgap(l-K) > ^| gap (I-K). 
We have to study the eigenfunctions / of P such that the corresponding <I> j satisfies 

*/ = 0. (3.17) 

Then Theorem 2.1 implies that $/ = ^4n for some AeM. On the other hand every tp£ 
can be decomposed as 

<pe = atn + (pt 

with ci£ = (<pe,n) and ((pi,n) = 0. Note that in view of (3.14) and the conservation law 
Yle^ivi) = 0) there is no restriction in assuming J2e a e = 0- Therefore 

A (n, n) = n) = (n, n) at = 

Thus if solves (3.17) then necessarily <3?j = 0. It remains to study the latter case in 
detail. 

Typical examples for which <£j = are obtained by choosing each tpg proportional to 
n. Call A the class of all / £ L 2 (v) of the form (3.14) with ipg = aifi with arbitrary 
a := {ai, . . . , ol} G C l . A simple computation gives 

(1-^)/=!^/, f€A (3.18) 

By (3.18) we may then restrict to the orthogonal complement A ± to prove our claim 
(3.13). We are going to prove that there exists a finite Lq = Lo(q) such that if / € A 1 - 
and = 

(1 - A) > ^| , VL > L . (3.19) 

Let us first check that if / € A 1 - then the corresponding tpg are all orthogonal (in H) to 
the number of particles: 

<^,n)=0, £=1,...,L (3.20) 



XL -I 



L 



K + 



K + 



L-l 
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Indeed / € A 1 - means 

= Y a k is((pi(r]i)n k ) = Y a k[Y^ e ' Kn) + ((p k , n) J 

L 



k,i 
L 



~T7 ak [ ~ n ) + L (Vk, n) = j-^j Y a k (fk,n) 



Here we are using (2.1) and <3?j = 0. The above identity implies (3.20) in view of the 
arbitrariness of a. Using again $j = we compute 



K/ 2 ) = ^2"(fk(Vk)feM) = Y [52(<Pk,K<Pt) + {<Pk,<Pk) 

k,t k £^k 

= Y [ {<Pk, K$ f ) + {<p k , (1 - K)<p k ) ] = Y ifk, (I " K)tp k ) 



A similar computation yields 

K/(i - P)f) = \Y ( J - K ^ L - V 1 + 

k 

Writing <pt = (I — K) 1 / 2 ^^ one has 

(i - a) y (<pt, <pe) = \ Y ^ L - 1 ) 1 + K )&) ■ 

i i 

Now observe that (<pe,n) = 0. This follows from (3.20), (2.1) and the self adjointness of 
K. By Theorem 2.1 we then infer 

for some uniform constants S > 0, k < oo. We conclude that 

L - 1 



(1-A)^ 



k L-2 



L L 2 + s " L - 1 
for L large. This finishes the proof of (3.13). 



□ 



3.2. A remark on Bernoulli Laplace model of diffusion. We observe that the strat- 
egy of the above proof may be used to compute in a simple way the spectral gap for the 
so-called Bernoulli-Laplace process, see [13], [12] and references therein. The latter can 
be seen as an exclusion process on a complete graph: there are L sites with exchanges 
allowed between any couple of sites i,j = 1,...,L and with uniform rates. The 

configuration space is Q ■= {0, 1} L and the measure v is the product Bernoulli measure 
on Slo conditioned to the event ^ a% = N . The Dirichlet form is then defined by (3.4) 
with Dij replaced by 



1 



Eij(f)=2"[(Vijf) 2 ] 



a k k^i,j 
a>i k = j 
yctj k = i 



(3.21) 
(3.22) 
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We are going to show that j(L) = i, where 



7(L)= sup sup ^1 (3.23) 
As in the proof of Proposition 3.1 - see (3.12) - we obtain 

- P)f) < - 1 W, /) • (3-24) 

where again P is defined by (3.10) with Th the cr-algebra generated by au € {0, 1}. The 
analysis of the spectrum of P is much simpler now. Indeed any / of the form (3.14), with 
v(f) = 0, here must be of type 

. N 

/(«) = 2^ aeae ' a e :=a e - — . (3.25) 

As in (3.18) a simple computation shows that if / is given by (3.25) then (I— P)f = jjz^f- 
It follows that any / £ L 2 (v) such that = satisfies 

„(/(I-P)/)>£z2„(/2) (3.26) 

By (3.24) we then obtain 

7(L) < 7(^-1) < ••• <7(2) ^^3. (3.27) 
On the other hand simple computations show that whenever N = 1 one has 

£(/,/) = 2Var v (/) 

for any / G L 2 {v). In particular 7(2) = 1/2 and 7(L) ^ 1/2 for all L ^ 3. Together 
with (3.27) this shows that 7(L) = 1/2. A more detailed spectral analysis of this model 
including all the eigenvalues was obtained in [13] by a different technique. 
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4. Quantum XXZ Hamiltonian 

Given S € \N, H G N, consider the Hilbert space ft = <8>^ =1 C 25+1 . The spin-S XXZ 
chain on [1, H] n Z with fcinA; boundary conditions is defined by the operator 

w (s) = Ewi3 + i. 
= 52 - A_1 (^^+1 + - ^^+1 + s ^ - A ~ 2 (^+1 - ^) • 

Here S l h , i = 1,2,3, are the spin-S operators (the 25 + 1-dimensional irreducible rep- 
resentation of SU(2)) at every h, and the constant S 2 has been added in order to have 
zero ground state energy. The parameter A £ (l,oo) measures the anisotropy along 
the third axis. The kink boundary condition is obtained through the telescopic sum 
Sjj — Sf = J2h=i {^h+i ~ ^h) an d the pre-factor Sy/l — A~ 2 is chosen in order to obtain 
non-trivial ground states describing quantum domain walls (see [1], [14] and references 
therein). 

We choose the basis of $j labelled by the 25+1 states of the third component of the spin 
at each site and we write it in terms of configurations 

m = (mi, . . . ,m H ) € {S, -S + 1, . . . , S - 1, S} H =: Q s 

so that \m) = ®^ =1 \rrih) stands for the generic basis vector in fj. With these notations, 
and introducing the stair-operators S ± = S 1 ±iS 2 , the action of S l , i = 1,2, 3, is given by 

S 3 h \m h ) = m h \m h ) , S± \m h ) = c±(S, m h ) \m h ± 1) . (4.2) 

c±(S,m h ) := y/(STm h )(S±m h + l) 

The action of TiS s ^ is explicated by rewriting the pair- interaction terms as 

«iS+i = S 2 ~ ( 2A ) _1 ( S hS h+1 + S h S+ +1 ) - S 3 S 3 h+1 + 5V1-A-2 {S 3 h+1 - S 3 ) (4.3) 

4.1. The spectral gap. The Hamiltonian commutes with the total third component 
of the spin 

H 

&tot = • 
h=l 

We shall divide the space S) into sectors fi n , n G {—SH, —SH + 1, . . . , SH — 1, given 
by the eigenspaces corresponding to the eigenvalue n of S 3 ot . It is known [1] that for each 
n there is a unique (up to multiplicative constants) vector ip n £ fj„ such that H^tpn = 0, 
which is given by 

meQs- 
Y, h m h =n 

Here g£ (0, 1) is the anisotropy parameter linked to A by the equation 

A=i( 5 + g- 1 ). (4.5) 
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The ground states ip n are interpreted as describing an interface profile, [1, 11]. A funda- 
mental question associated to the stability of such "quantum interfaces" is the positivity 
of the spectral gap [9, 10]. The latter, denoted gap(H^)), is defined as the energy of 
the first excited state, i.e. the first non-zero eigenvalue of the non-negative operator 
Recently this question was studied in great detail in the paper [10] by both analytical 
and numerical means. One of the main results of [10] is a proof of the fact that for every 
S G |N, gap(H^) is positive uniformly in H. Furthermore it was conjectured on the 
basis of numerical analysis that gap(TC^) should grow linearly with S. We prove the 
following bounds. 

Theorem 4.1. For every A G (l,oo), there exists 5 > such that 

5S ^ gap(H {5) ) ^ S^S 

for all S G and all H ^ 2. 

In order to prove Theorem 4.1 we shall establish the following unitary equivalence. Let 
L = 2S and N = SH + n, recall the definition of L 2 (S7, i>), of the subspace of horizontally 
symmetric functions and of the generator C introduced in (1.12). The measure v in (1.11) 
can be written, using (4.4) with m = to — S, as: 

v{u) = 1[^(uj-S)\ 2 Z= Y, [V^-S)] 2 - (4.6) 

Y,h UJ h=SH+n 



For any ip G $j n we also write 



<P 



m€Q S - 
T, h ™ h =n 



Finally we set <p(u) = ip(oj — S). Then the transformation 

l f( m ) -> rtr^: =: \Pnf] M > uj = m + S 

maps unitarily S) n into L 2 (f2,z>). 

Lemma 4.2. For every n G {-SH, -SH + 1, . . . , SH - 1, SH} 

U n HW<p = ~CU n <p, V G ?) n (4.7) 

Proof. From (4.3) we compute 

ft^h+i) l m ) = ( S2 ~ m h m h+i + S\A - A- 2 (m h+1 - m h )) \m) 

- (2A)- 1 c + (S,m h )c-(S,m h+1 )\m + > h ) - (2A)- 1 c^S,m h )c + (S,m h+1 )\m^ h ) . 
Here we are using the notation 

mh' h! h, h + 1 

m h' h = ' m h ^ 1 h' = h 

m h+1 =F 1 h' = h + 1 
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Therefore 

(2A) [H { ^ h+1) ip] (m) = T(m h , m h+1 ) (p(m) 

- c+(S,m h )c-(S, m h+1 )ip(m + ' h ) - c-(S, m h )c+(S, m h+1 )ip(m^' h ) 

with 

T(m h ,m h+1 ) = (2A)[S 2 - m h m h+1 + Sy/l- A~ 2 (m h+ i - m h )] . 
Now a computation shows that 

T(m h , m h+ i) = w +th (co) + w- th (u) , to = S + m , 

with w±f l the rates defined in (1.12). Another computation shows that 



c + (S,uj h - S)c-(S,u; h+1 - S) J v(yU) - = w +!h (u) 

V v \ u ) 



and similarly 



c-(S,u h - S)c+(S,u h+1 - S)J U<Kl f - = . 

V 



We have then obtained 

(2A) [U n H%] h+l) v] (u,) = {w +A (u) + w_, h {u)) [U n <p] {u) 

- w + , h (oj) [U n <p] (uj + ' h ) - w- h {u) [U n <p] (u~' h ) 
and the lemma follows. □ 

We are now able to finish the proof of Theorem 4.1. Recall (1.7). Since here L = 2S, we 
readily infer from Lemma 4.2 the estimate 

where j(2S,H) is defined by (1.8). The bound 7 ^ 5~ l is the content of Theorem 1.1. 
On the other hand to prove gap^^) ^ <5 _1 S we may use the following simple argument 
which says that in each sector Sj n there are excited states with energy bounded by 5 -1 5. 
Choose f(oj) = 0Jh with ho = [p], p = N/(2S). A simple estimate shows that 

£(/(-£)/) < q^u{uj ho + u ho+1 ) +qv(u ho + w ho _i) ^ k S {I A p) 

for some finite A; = k(q) < 00. On the other hand using the estimates in [7] it is possible 
to check that Varp(/) ^ 5 S (1 A p) for some S = 5(q) > 0. By the variational principle 

Kap(W (*))< * *(/(-*)/) <fc g 
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5. Energy gap above the diagonal interfaces of the XXZ model 

In this section we study a higher dimensional quantum XXZ Hamiltonian which is 
sometimes used to model a tilted interface ([11, 4, 7]). In order to avoid complicated 
notation we shall work in a two-dimensional setting. We later observe that our results 
actually hold without modifications in any dimension. 

We begin by defining a cylindrical region (in two dimensions simply a rectangle) with axis 
along the diagonal. Given two integers R and H we define 

r = T RyH = {x G 1? : -R < xi - x 2 < R , 1 < xi + x 2 < H} (5.1) 

We write l x = x\ + x 2 for the distance of a site x form the line X\ = —x 2 . A bond is an 
oriented pair b = (x,y), with x,y G Z 2 such that |xi — y±\ + |x2 — 2/2 1 = 1 (in particular 
£y = i 1). We call B = Br^h the set of bonds b = (x, y) with x, y G T and i y = £ x + 1. 
For any S 1 G the anisotropic spin-5 Hamiltonian in the region T with kink boundary 
conditions is defined by 

wf> = S 2 - A" 1 (S^J + S 2 x Sl) - SlS} + SVl ~ A- 2 (S, 3 - S 3 X ) , b=(x, y) 

with S x , k = 1,2,3, the spin-,!? operators at site x. This is the higher-spin analog of 
the spin-^ cylindrical models considered in [4, 7]. As usual we consider the Hilbert space 
£)r = <8'a;erC 2 ' s+1 with the basis labelled by configurations 

\m) = ® xe r \m x ) m G {-S, S} r =: Q r ,s ■ 

Clearly the total third component Sf ot = Ylxer ^x ^ s conserved and we may divide f)r into 
sectors $)r,m n G {— S\T\, . . . , S\T\} according to the eigenvalues of Sf ot . Following [1] we 
know that in each such sector there is a unique ground state tpr,n given by 

Vr,n = ^2 ipr(m) \m) 

meQr,s- 

*•<■»> -n^GfJ- < 5 - 3 > 

As before q = q(A) is defined by (4.5). The above ground states have zero energy and we 
call gap(W^p) the first nonzero eigenvalue of T~l^\ 

Our main result here is a generalization to higher-spin models of a theorem we proved in 
the case S = \ , [7] . 

Theorem 5.1. For every A G (l,oo), there exists 5 > such that for all S G ^N, all 
R, H G N 

5 S R- 2 < gap(wjf } ) < (T 1 5 iT 2 . (5.4) 

Remark 5.2. It will be clear from the proof that exactly the same estimates hold (except 
that now the constant 5 in (5.4) may depend on the dimension) in the case where T is 
replaced by a cylinder in d + 1 dimensions with axis along the 11 • • • 1 direction and basis 
given by a d-dimensional hypercube of side R. 
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Proof. As a first step we establish a unitary equivalence in the spirit of Lemma 4.2. Namely 
let Clr = (25 + l) r denote the set of configurations w x G {0, . . . , 25}, x G T. Following 
our previous analysis we may interpret oj x as the number of particles at site x. Given 
N G {0,1,... ,25|r|} define " 

&r,N = G &r ■ ^ Ux = N } 

Consider the probability measure on Or,iV given by 

*rwM = -^[^r(^-5)] 2 , Zryv = V [* r (^-5)] 2 . (5.5) 

Then as in Lemma 4.2 we obtain that 

f( m ) 7= 1 7 ; ^(^ - 5 ) = : [ u r,n<f] M , cj = m + 5 
V^T,at(w) 

maps unitarily each sector f)r,n into L 2 (£Iy,n, £r,7v)> N = S\T\ + n. Repeating the com- 
putation leading to (4.7), for every bond b G B we have 

U T , n Wj 5) f = ~JK i< V tf + n"V 6 -/} , / G L 2 (f2r,iv, *>r,iv) (5.6) 
with the notation 

u>b(v) = q' 1 uj y {2S - uj x ) , w^iuj) = quJ x (2S - u y ) , b= (x,y) , 

V±/H = /(^' ± )-/M, <4*:=lu; x ±l z = x 
Let us define the Markov generator Gr,s by 

0h,s/M = £ + nHV 6 -/H} . (5.7) 

beB 

Then is symmetric in L 2 (&r,N, vt,n) an d by (5.6) we have the unitary equivalence 

-Hf^-^Qw (5.8) 

Set r2 r (s) = {0, l} r(S) , where is the 3-dimensional region 

r( s ) = {(i,x) : i = 1,... ,25; x G T) . 

Define then f2 r (s) N as the set of a G O r (s) such that Ylu x)gr( s ) a (i,:r) = AT. Consider the 
probability measure vn on fi r (s) tv defined by 

W Mer(f) a en r(S) jv (j,y) G r(s) 

Consider the subspace 5r C L 2 (Q r (s) n ,vn) of symmetric functions defined by: f(a) = 
f(a w,x ) for all a, all x G T and all permutations ir G T^Sj with a 71 "^ defined as in (1.9) 
replacing h by x. As in section 1, vy,n can be looked at as the marginal of vn on the sums 

2S 

w x = £ a (i,x) > a ^ ^F(S),N 

1=1 
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and Sr is identified with L 2 (£7r,7v, vt,n)- Then Gr,s niay be identified with the restriction 
to Sr of (2S)Gr,s, the Markov generator defined by 



25 25 



Gn,sf{a) = j^YY Y c (<,*);(j.w)( a ) v (<.*);C?.v)/( a ) ( 5 - 10 ) 

i=l j=l ( X ,i/)€B 

c {i,x);{j,y) ( a ) = 9 a < i -)- a W-») , V ( ^ );0 - = /(a^W*)) - /(a) . 



Here as usual q^^)^^) denotes the configuration after the exchange between (i, x) and 
(j, y). At this point we have obtained a unitary equivalence 



(5.11) 



when the right hand side above is restricted to Sr- Notice the analogy of (5.10) with 
the process introduced in section 1, see (1.4). However there is an important difference 
(which will be seen in a moment to be responsible for the R~ 2 factors in (5.4)): while 
particles diffuse asymmetrically in the 11 direction just as it happens for (1.4) along the 
vertical direction, we have in (5.10) in addition an essentially symmetric diffusion along 
the orthogonal direction (given by the lines i x = constant). We are now going to take care 
of these facts. 



Observe that Qr^s is symmetric in L 2 (Q r (s), vn) with Dirichlet form 

25 25 

s R ,s(f) = VN{f(-g R ,s)f) = ^ EE%>(/) 



(5-12) 



i=l j=l 



E (i,j)(f) '■= n Y VN \ c (i,x);U,y)( v (hx);U,y)fY 



(x,y)eB 

For any x 6 T let t x = x\ — X2, so that x is uniquely determined by a couple (t, £), 
—R ^t^R, l^l^H. We may then write 

1 21 

E (iJ)(f) = 2 Y x(t,t)"N [C(i,tfl-,(j,8j+i)(V(i,tfl-,(j,8j+i)f) \ (5-13) 

t,se[-R,...,R]: 1=1 
\t-s\=l 

where x(i,^) = 1 if there exists x G T such that x\ — X2 = t and x\ + x 2 = t and 
x(t, t) = otherwise (notice that if x(^) = 1 then xi s ^ + 1) = 1 for s G [-R, ■■■ ,R] 
with \s — t\ = 1). Using the bounds q ^ c ^ q^ 1 on the rates and the properties of the 
measure un it is not difficult (see Lemma 2.6 in [7] for similar computations) to show that 
there exists k = k(q) < oo such that given arbitrary t,s£ [-R, ■ ■ ■ , R] with t ^ s — 1 and 
£ G [l,H] one has 



x(t,l)x(s,£+l)vN [c(i,t,ey,(j,s,e+i)(V(i,t,ey,(j,s,e+i)fy 



s-1 



r=t 



(5-14) 
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c(i,t,ey,(j,s,e+i) (V(i,t,ey,ti,s,e+i)fY 



From (5.13) and (5.14) we easily obtain 

H-l 
t,se[-R,...,R] t=\ 

< kR 3 E (hj) (f) (5.15) 
Recalling (5.12) we may summarize the above estimate with the following statement 

£ R ,s(f) > SR~ 2 V R , s (f), (5-16) 
1 



2SR 



E 2 D (i,t);(j,s)(f) , 



i,je[i,...,2S] t,se[-R,...,R] 



H-l 



We are now able to conclude thanks to Theorem 1.1. Indeed, the Dirichlet form V^s 
defined above is a special case of the one appearing in the Theorem, when we set L = 2SR 
if R is even and L = 2S(R + 1) if R is odd. We then have that 

V R , s (f) > <5Var^(/) 

for some uniform constant 5. From (5.16) and (5.11) we see that the energy of excited 

states in each sector lies above 5 S R~ 2 for some uniform constant S > 0. The desired 

( s) 

lower bound on gap (H%') follows at once. The reverse estimate is much easier and it can 
be obtained again from (5.11) by modifying slightly the reasoning in Proposition 6.1 of 
[7]- 

□ 
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